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Driving a double-quantum-well excitonic intersubband resonance 
with a terahertz (THz) electric field of frequency lot Hz generated ter- 
ahertz optical sidebands lo = lothz + unir on a weak NIR probe. 
At high THz intensities, the intersubband dipole energy which cou- 
pled two excitons was comparable to the THz photon energy. In this 
strong-field regime the sideband intensity displayed a non-monotonic 
dependence on the THz field strength. The oscillating refractive in- 
dex which gives rise to the sidebands may be understood by the for- 
mation of Floquet states, which oscillate with the same periodicity 
as the driving THz field. 
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Electromagnetically-induced coherent quantum effects in 
semiconductor systems include Rabi oscillations in donors [|lj| 
and quantum dots [|[, and electromagnetically-induced trans- 
parency in quantum wells [j^]. In these effects, to induce co- 
herence the driving field must overcome dephasing processes. 
Thus, the Rabi frequency fiE must be greater than the dephas- 
ing rates 7: [iE > hj, where fi is the dipole moment and E 
is the electric field strength. 

If the field strength is increased still further, a strong-field 
regime is encountered as the Rabi energy becomes compara- 
ble to the photon energy: /j,E > hw. Previous studies into this 
regime have looked at atoms in microwave cavities and 
strong laser fields [|J|. Since atomic excited states are near the 
continuum, strong-field effects manifest themselves as non- 
monotonic multi-photon ionization rates as a function of field 
strength. 

However, the atoms are ionized, meaning the bound states 
do not exist after the microwave or laser field is turned on. 
Multi-photon ionization is an inescapable result of driving ex- 
cited atomic levels which get closer and closer together. 

There is a rich body of theoretical predictions on the effect 
of strong fields on the bound states of a quantum system l^js]]. 
In quantum wells, upper states get further and further apart. 
Thus the strong-field condition can be satisfied by resonantly 
driving the lowest subbands of a quantum well at THz fre- 
quencies. Yet the quantum well is far deeper than the Rabi 
energy, so the bound states still survive. 

We describe experiments in which a near-infrared (NIR) 
probe laser beam is mixed with a terahertz (THz) pump beam 
in a gated, asymmetric double quantum well (DQW). The 
THz field couples to an excitonic intersubband excitation 
while the NIR field couples to an excitonic interband exci- 



tation. Applying a DC voltage to the gates brings the inter- 
subband transition into resonance with the THz field, and the 
interband transition into resonance with the NIR field. When 
these resonance conditions are met the NIR probe is mod- 
ulated resulting in the emission of optical sidebands which 
appear at frequencies Lo S ideband = u NIR + nuj THz where 
unir (uthz) is me frequency of the NIR (THz) beam and 
n = ±1, 2, . . .. Strong-field effects manifest themselves as 
non-monotonic sideband intensities as a function of THz field 
strength. 

The sample consisted of an active region, gates, and a dis- 
tributed Bragg reflector (DBR). The active region consisted 5 
periods of DQW, each consisting of a 120 A GaAs QW and 
a 100 A GaAs QW separated by a 25 A Al0.2Ga0.sAs tunnel 
barrier. Each period was separated by a 300 A Alo.3Gao.7As 
barrier. The active region was sandwiched between two gates, 
each consisting of a Si delta-doped 70 A QW with carrier den- 
sity Ri 1 x 10 12 cm' 2 . The gates were separated from the 
active region by 3000 A Alo.3Gao.7As barriers. Since the gate 
QWs are much narrower than the active DQWs, the gate QWs 
are transparent to both the THz and NIR beams. The gated 
DQWs were grown on top of a DBR which consisted of 15 
periods of 689 A AlAs and 606 A Alo.3Gao.7As. It had a low- 
temperature passband nearly centered on the low-temperature 
bandgap of the DQW, making it about 95% reflective for the 
NIR probe beam and sidebands. 

We etched a mesa and annealed NiGeAu ohmic contacts to 
the gate QWs. The sample was then cleaved into a 1 mm wide 
strip, 8 mm long. We then cleaved a 400 /im-thick wafer of 
crystal sapphire substrate material into a 1 x 5 mm strip, ori- 
ented with the optically slow axis along the long dimension. 
The sapphire was mechanically pressed against the surface of 
the sample with a beryllium copper clip as shown in Fig. [jj 
At cryogenic temperatures sapphire is index-matched to GaAs 
and transparent at THz wavelengths. At the same time, sap- 
phire is transparent to the NIR probe. Pressing the sapphire 
against the sample forms a rectangular dielectric waveguide 
with half of the waveguide defined by the sample substrate and 
the other half defined by the sapphire. The epilayer containing 
the active region lies in the middle of the heavily over-moded 
waveguide. This results in much higher fields at the active re- 
gion than would be possible without the sapphire, where the 
active region would be at the edge of a dielectric waveguide. 

The sample was cooled to 21 K in a closed-cycle He cryo- 
stat. The 1 kW THz beam from a free-electron laser is polar- 
ized in the growth direction, propagated in the QW plane and 
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focused by a 90° off-axis parabolic mirror (F/l) into the edge 
of the waveguide. The maximum THz electric field strength at 
the focus is estimated to be between w 5 — 20 kV/cm, where 
our uncertainty comes from uncertainty of the exact size of the 
mode in the waveguide. The THz power can be continuously 
varied by a pair of wire-grid polarizers. 

NIR light from a continuous-wave Tksapphire laser is 
chopped by an acousto-optic modulator into 25 /is pulses 
which overlap the 5 /is FEL pulses at maximum repetition 
rate of 1 .5 Hz. The vertically -polarized NIR beam propagates 
normal to the THz beam, and was focused (F/10) at a power 
density of » 5 W/cm 2 to the same small interaction volume 
in the sample. The reflected beam, sidebands, and photolumi- 
nescence (PL) are analyzed by a second polarizer, dispersed 
by a 0.85 m double-monochromator, and detected by a photo- 
multiplier tube. 

Previous work [^,|l0|] showed that sideband generation is 
enhanced when the NIR field resonantly couples the vacuum 
state |0) with an exciton |1), and the THz field resonantly cou- 
ples with the intersubband transition between excitons 1 1) and 
|2). The exciton states can be labelled as E^HH^X, indicat- 
ing an exciton consisting of an electron in conduction sub- 
band /t and a heavy hole in valence subband v. Changing the 
DC electric field tunes the excitonic intersubband transition 
into resonance with the THz field, as shown in Fig. |2| The 
sideband resonance we focused on is E2HH2X-E2HH3X, a 
peak assignment made by comparing low-field sideband spec- 
troscopy with a nonlinear susceptibility calculation for exci- 
tons Jio| . 

FigTTfl shows the dependence on THz field strength of side- 
band generation at various THz frequencies. Each curve was 
taken at a gate bias near the E2HH2X-E2HH3X resonance 
JTo|]. The most striking feature is the non-monotonic behavior 
in the strong-field regime when the Rabi energy is comparable 
to the photon energy. Clearly at lower frequencies the power 
dependence rolls over at a lower field strength. 

By sitting at the peak NIR frequency for E2HH2X- 
E2HH3X and varying the DC gate voltage we took THz power 
dependence scans at various THz detunings E2—E1 — Tiujthz 
where E2 (Ei) is the energy of the upper (lower) exciton state. 
The results are shown in Fig. [|. As the detuning is varied, so 
does the shape of the THz field dependence. 

The power dependence cannot be explained by a nonlinear 
susceptibility, which is inherently a low-field theory because 
it relies on a Taylor expansion about the field strength JTT|]. 
Saturation effects [ |l2] ] in a nonlinear susceptibility can only 
come about when contributions from virtual transitions ini- 
tiated from excited states |1) or |2) are comparable to those 
initiated from the vacuum |0). However the populations of ex- 
cited states |1) or |2) are never significant compared with |0) 
in our undoped sample and weak NIR beam. Therefore a non- 
linear susceptibility can only predict a linear (quadratic) de- 
pendence of sideband intensity on THz power (field strength). 

The following simple, phenomenological, three-state 
model captures qualitative features of the dependence of side- 
bands on THz power. The Hamiltonian Ho has eigenstates 










°\ ( 










Zl2 \ x = 




Vo 


Z21 


Z22 J \ 



(j>o(x,z), cf)i(x,z), and (f>2(x,z) with eigenenergies E = 0, 
Ei, and E2, respectively. Near the experimental resonance, 
the effect of the DC voltage is to tune $2 while the other <f>i 
remains relatively unchanged [|l0|], The two exciton states 
4>i{x, z) and <f>2{x,z) are coupled by the z-dipole operator 
for the DC and THz fields. The two upper states <fri (x, z) and 
<j>2(x, z) are coupled to the ground state by the x-dipole oper- 
ator for the NIR field. 

Therefore, the x and z matrices are 

xoi X02 
x w 
X20 

All the nonzero terms x a p = (<j> a \x\(f)p) and z a p = 
(<j> a \z\<j)p) are set equal to unity. The Hamiltonian is 

H = Ho + zE^ cos u>t + xXEq cos fit 

where E u represents the strong THz electric field strength 
and Eq represents the weak NIR electric field. We solved 
the Hamiltonian nonperturbatively within a Floquet formal- 
ism for A = 0, while the weak probe was added later using 
time-dependent perturbation theory. 

The solutions to the Schrodinger equation for the time- 
periodic Hamiltonian have the form Jl3| ] 

<fi(z, t) — e~ l€it / h iii(z, t), u(z, t) = u(z, t H ) 

OJ 

where Ui(z, t) has the same periodicity as the driving fre- 
quency and and e; is called the quasienergy . The states 
ifi (z, t) are called Floquet states and are mathematically anal- 
ogous to Bloch states for a spatially-periodic Hamiltonian. 
Meanwhile the quasienergies are analogous to the crystal mo- 
menta of Bloch theory. 

The spatial-dependence of the wavefunction <pi(z, t) can be 
expanded in terms of the original eigenstates <p a (z), and the 
time-dependence can be Fourier-expanded 
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The coefficients c l a n are what will determine the part of the 
index of refraction which will oscillate at particular multi- 
ples of oj. Solving for these coefficients is the key to un- 
derstanding the power dependence of the sideband. These 
Floquet coefficients are closely related to the photon-assisted 
tunnelling sidebands which appear in the irradiated current- 
voltage curves in superconducting weak-link junctions Jl4| ] 
and coupled quantum wells jl5|]. 

The time-dependent Schrodinger equation can be cast in 
terms of the Floquet operator F = H ' 
written as a matrix of the form [f§] 
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where [Ho] and [/i] are matrix representations of the opera- 
tors Hq and /i in the \<j> a ) basis. In practice the Floquet matrix 
must be truncated up to ±N photons. N can be made arbitrar- 
ily large for an arbitrarily precise result at high field strengths. 
The solution to the Schrodinger equation for a time-periodic 
Hamiltonian reduces to finding the eigenvalues and eigenvec- 
tors of the Floquet matrix Other than the ±N truncation 
the solution involves no perturbation or rotating-wave approx- 
imations. Here we used N — 16 photons. 

The Floquet solutions are labeled (pi(x, z,t) with i = 
0,1,2, and expanded as in (jl|). Their expansions are explicitly 
expressed for clarity: 

(f (x,z,t) = <j) Q (x,z) 

Mx, z, t) = e*i*/» J2 e- inut [cl,n4>i& *) + c^M^ z)] 

n 

y 2 {x, z, t) = e^/ h e~ m ^[cl n M^ z) + <% in <hfa *)] 

n 

The ground state <j>o(x, z) is not coupled by the the strong 
field to the 2-d subspace spanned by \<pi) and [4>2). Therefore 
the coefficients for the Floquet state cpo{x, z i t) vanish except 
c[] q = 1. In other words, the ground vacuum state remains 
untouched by the intense THz field. 

The dipole response of the Floquet states to the weak NIR 
probe is calculated perturbatively. The state of the system with 
both driving fields ip(x, z, t) can be expanded in terms of the 
Floquet states <fi(x, z,t) to first order. Discarding second or- 
der terms and anti-resonant contributions, the dipole expecta- 
tion is given by 



iQI 



E 



((po\x\ipi) ((pi\x\ip ) 
e, — hQ 



(3) 



The form of (||) is the same as the linear susceptibility of an 
undriven system, except the states ipi are oscillating Floquet 
states instead of stationary states. Explicitly expanding the 
Floquet states in the numerator gives us an expression for the 
polarization 
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The +1uj sideband is given by the Fourier component of 
x(t) oscillating at the frequency fl+lui. Also, for the resonant 
condition illustrated in Fig. |[ Tiui rs ei, so we keep only the 
i = 1 term in the sum (^). This condition is satisfied when 
n — m = 1. Thus we obtain an expression for the sideband 
polarization 



X sideband (^) 



€l — TiQ 



E c l+i,o ! 4i,i3 x opx a o (5) 



The intensity of the sideband is proportional to x 2 sideband . 
This is plotted vs. field strength for various detunings in Fig. 



|[ The detuning parameter d is the level spacing normalized 
by the photon energy of the strong field, d — ^j^f 1 ■ 

Our model captures the rollover of the resonant power de- 
pendence up to field strengths of around w 1.5. Given a 
calculated excitonic intersubband dipole of ^ = 12 nm JlO|], 
the strong-field condition is met at a THz field of 5 kV/cm. 
However, experimentally the sideband never completely dis- 
appears, a striking prediction of the theory. It is unlikely that a 
3-state simplification is entirely valid because there is a spec- 
trum of other exciton states as well as an electron-hole contin- 
uum that may provide significant off -resonant contributions. 

To summarize, we drove an quantum-well excitonic inter- 
subband transition with an intense THz laser field in a regime 
where the Rabi energy was comparable to the photon energy. 
The nonlinear mixing between the THz pump and a weak NIR 
probe displayed a nonmonotonic dependence on THz power, 
which could not be explained by a conventional nonlinear sus- 
ceptibility. Rather, the excitons were dressed by the THz field, 
a process which we described in a simple three-level model in 
which the interaction between excitons and the THz field is 
solved non-perturbatively wit hin the Floquet fo rmalism. 

This research is funded by frisf-dmr 0070083 . 
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FIG. 1. Processed waveguide device and optical coupling scheme. 
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FIG. 2. The dc electric field dependence of sideband at 
Tiwnir = 1547 meV and luthz = 1.5 THz (6.2 meV). The inset 
illustrates the resonant condition in which the NIR field is resonant 
with the E2HH2 exciton and the THz field resonantly couples the 
E2HH2 and E2HH3 excitons. 
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FIG. 4. THz field strength dependence of n = 1 sideband for 
o!jv/i?=1548 meV, ujthz = 1.5 THz (6.2 meV) at various DC elec- 
tric field detunings from the resonance condition illustrated in Fig. 
^. The absolute THz electric field scale is accurate only to within a 
factor of 2. 
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FIG. 3. Resonant THz field strength dependence of sideband gen- 
eration at 1.5 THz, 2.0 THz, and 2.5 THz (6.2 meV, 8.2 meV, 10.4 
meV). The absolute THz electric field scale is accurate only to within 
a factor of 2. 




Field strength irE/Tico 
FIG. 5. Sidebands calculated by evaluating the square of (Q) at 
various detuning parameters LEFT: d — E2 r ~ u f 1 = 0.8, 1.0, 1.2, 
RIGHT: d = E2 ~ E1 = 0.2, 0.4, 0.6. 
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